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LIST OF SYMBOLS 



s. 

fi 


Symbols 

b 

c 

S 

M 

N 

x,y,z 


Greek 

a 

Subscripts 

cp 

i 

J 

k 

I 


Description 

Wing span 
Local chord 
Root chord 

Number of spanwise strips plus one 
Number of bound elements 
Wing rectangular coordinate system with x 
in the streamwise direction and y to the 
right 

Angle of attack deg 

Control point 

Chordwise bound element number 
Spanwise strip number 
Chordwise bound element number 
Leading-edge 


Dimensions 
m (ft) 
m (ft) 
m (ft) 

m (ft) 






1. IHTROKJCTIOH 

'Hiis docioent describes in detail ths necessary information for using a 
computer program to predict distributed and total aerodynamic characteristics 
for low aispect-ratio wings with partial leading-edge separation. This program 
is based on the numerical method developed in reference 1. The flow is assumed 
to be steady and Inviscid. The wing boiindary condition is formulated by the 
Quasi-Vortex-Lattice method. The leading-edge separated vortices are 
represented by discrete free vortex elements which are aligned with the local 
velocity vector at mid-points to satisfy the force free condition. The wake 
behind the trailing-edge is also force free. The flow tangency boundary 
condition is satisfied on the wing, including the leading- and trailing-edges. 

The program is restricted to delta wings with zero thickness and no 
camber. It is written in Fortran language and runs on CDC 6600 Computer. 
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COIPUTER PROGRAM DESCRIPTIOll 


2.1 PROBLB! D^IHETIOH 

In steady symaetric flight at a high angle of attack, the flow over a thin 
low aspect-ratio highly sweptback wing separates along the leading-edge and the 
tips. In the following, only delta wings are considered. The wing can be 
represented by a bound vortex sheet, across which there exists a pressure 
difference, and the separated flow along leading-edges by force free vortex 
sheets, across which there is no pressure difference. In the present method, 
the Quasi-Vortex-Lattice method (reference 2) is used to simplify the induced 
velocity expressions due to the bound vortex sheet and discrete force free 
vortex elements for separated vortex sheets. The following boundary conditions 
sure imposed on the flow model: 

a. The flow must be tangential to the wing surface. 

b. The leading-edge boundary condition and the trailing-edge Kutta 
condition are to be satisfied. 

c. The vortex elements over the wing and wake behind the trailing-edge 
are force free. 

This is a non-linear problem because the strengths of the wing bound vortices 
and free vortices, and the locations of the free vortex elements are unknown. 
Thus, the problan is solved by an iterative method. 

2.2 PROGRAM CAPABILITIES 

This computer program provides a theoretical method for determining the 
aerodynamic characteristics of low aspect-ratio thin delta wings without camber, 
with partial leading-edge separation. The following is a list of the aerodynamic 
characteristics the program calculates: 

a. Spanwise and chordwise ACp distributions 
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b. Spanvlse distribution of sectional lift, induced-drag and pitching 
mtanent coefficients. 

c. Total lift, induced-drag, pitching-nooent and leading-edge thrust 
coefficients. 

2.3 GEOMETRY DESCRIPTION 

The origin of the rectangular coordinate system is at the wing apex. 

The wing lies in the x-y plane and the x-axis is taken along the wing center-line 
The wing span is given by b and the surface area S. . 

2.3.1 WING GEOMETRY 

The location of bound- and trailing-vortex elements for a typical case 
are shown in figure 1. The x-location of bound elements is given by the 
cosine law and is illustrated in figure 1. 

Xf. = I (1 - Cos( ^ ~ ^ - tt)), (1) 

i = 1,2, N 

where is the leading-edge x-coordinate, c is the chord and N is 

the number of bound elfflnents in a chordwise direction. The spanwise 
location of trailing elements irf given by, 

Yj = ^ (1 - Cos(- ^ ^- g" tt) ), (2) 

j » 1,2, M 

where b is the span and M is the number of legs of trailing vorticity, 
which is one higher than the number of spanwise strips of bound elements. 

The locations of control points are given by. 
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( 3 ) 


k - 0,1.2, N 

^cp, * ¥ 

J 

J . 1,2, (M - 1) 


where 



“■» 'j 
respectively. 


are the leading-edge x-coordinate and chord at 


It has been found numerically that the aerodynamic characteristics 
depended on the nimiher of spanwise strips, l.e. M of equation (2). 
Therefore, a parametric study has been made to find a relation between the 
aspect ratio and the number of spanwise strips for reasonably accurate 
results (Fig. 2) (Section 3 of ref. l). It is to be noted that as the 
aspect ratio is decreased, the number of spanwise strips has to be increased. 
This is due to the fact that the spanwise variation of aerodynamic character- 
istics, such as pressure coefficient and thrust coefficient, is large for 
small aspect ratio wings. This study was performed by matching the lift 
coefficients obtained by using the present method to those obtained by using 
suction analogy (ref. 3) at one angle of attack. 


2.3.2 LEADING-EDGE VORTEX SYS'^M GE(»!ETRY 

The leading-edge vortex system is superimposed on the regular quasi-vortex- 
lattice grid. A typical vortex elanent is shown by points A through J in 
figure 3. These points are connected by a series of short straight segments. 

The initial location of these segments is shown by dashed lines and final 
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Number of spanwise strips, (M-l) 


Figure 2. Variation of nuaber of spanviae strips with aspect ratio 
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locution ^ 80114 linos. Those segments hove the follovlng characteristics; 

a. Points A ttnrough E lie aloi^ a ving trailing vortex elentent. 

Initially point A is one root ctord avay front the trailing-edge 
in the downstream direction and the line segments between A and 0 
are parallel to the axis of symnetry. The line segments between 
points A and B are of equal length. In the final converged position 
these 8e£Bents are aligned in the direction of the local velocity 
vector. The segments B-C and C-D are 0.1 long. B-C is allowed 
to move only in the vertical direction whereas C-D is fixed it the 
wing plane because the flow is tangential to the trailing-edge. 
Segment D-E is also fixed in the wing plane. 

b. Points E, F, G and H also lie in the ving plane. The location 
of segment E-F is ahead of the wing first bound element and 

is given by, 

£ 

*r ■ % * r - =‘”'2(rh7” '5.t) 

where the subscripts E and F refer to *he points under consideration. 
The above two equations are similar to equation (l). It is to be 
noted that sequent E-F its located at the first bound element for a 
grid of (K + l) bound elments in a chordwise direction. The 
se^ents F-G and C-H are of tlie same length and point G lies on the 
leading-edge. The segment G-H is fixed in the ving plane due to the 
leading-edge boundary condition. 
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c. nie initial location of point X it given 


Xj - Xy (6.a) 

yj ■ yp 

Zj ■ 0.1 Cp tan(22.5 - 0.5o) for o £ 15° (6.c) 

or z_ ■ 0.1 C_ tan a for o > 15° (6.d) 


where C_ is the root chord and a is the angle of attack, 
n 

Initially point J is one root chord .away frewj the trailing-edge. 

The segments between point I and J are of equal length and lie in 
a plane parallel to x-z plane. These segments are approxim&tely 
at a height of 0.'l*Cp above the wing plane. In the final conve 'ci u 
position all the sepnents between points H and J are aligned in the 
direction of the local velocity vector, 
d. The s«ni-infinlte sepnents from points A to infinity and J to 
infinity are straight and are parallel to the undisturbed 
free-stream direction. 

2.U SOUJTION PBOCEDURE 

The basic unknowns of the problem are the bound vortex density on 
the wing, and the strengths and the locations of the element.«5 of the 
leading»edge vortex system and the wake. The problan is nonlinear 
because the locations of the leading-edge vortex system and the wake are 
unknown a priori. Therefore, the problem is solved by the iterative 
process described below; 
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a. Prescribe the vortex lattiee for the ving surface, and the initial 
locations of the free elements over the ving and in the wake. 

h. By satisfying the ving hoiudary condition, -obtain the bound vortex 
density of the ving and the strengths of free elements. 

c. Calculate all the aerodynamic characteristics. 

d. Calculate the forces acting on the free elements over the ving 
surface. 

e. Adjust the free elanents of the leading-edge vortex system and 
the vake in the local velocity vector direction. 

f. Repeat steps b through e until a converged solution is obtained. 

The initial locations of the free vortex elements are assumed by 

letting them leave the leading-edge- in the undisturbed free-stream direction 
up-to a height of about ten percent of the root chord beyond which the 
elonents are parallel to the wing plane. Initially, all the elenents of 
the wake lie in the plane of the wing. In the iteration process, the 
force free condition is satisfied on the free elements from the root to the 
tip in t ,, down-stream direction. The eltments over the wing are adjusted 
before the elements of the wake. In the first iteration the segments over 
the wing are moved 100 percent according to the velocity computed at their 
mid-points. This movement is gradually reduced in steps of 90» 80 and 
75 percent in the next three iterations, after which it remains at 
75 percent (Section 2.5.2 of ref. l). The segments in the wake are moved 
only 50 percent in each iteration. Thus, exact force free condition is not 
enforced because whenever the free elements come close to each other they 
induce unreasonably large velocities because viscous effects are not 
included in the present theory. These large velocities increase the forces 
on the segments and induce fluctuations in their locations. 
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The solution is assumed to have converged if in tvo consecutive 
iterations the difference hetveen t)M total strengths of leading-edge free 
vortex elements is less than one percent and the absolute force acting on 
the free elen^nts is in the i^ighhorhood of a minimum^ Thus, an exact 
force free condition is not enforced as discussed in the previous paragraph. 
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3. IHPUT DATA FORMAf 

The following is the description of input data cards for this program. 
Card 1. Fomat (i6A5) 

TTL Any title identifying the case to he run. END in first 

three columns terminates the Job. 

Card 2. Format (6X3) 

NCW Number of chordwise lines (limited to nine) 

NSW Number of spamd.se lines (one higher than number of 

speuiwlse strips of pcmels, limited to twenty). It 
depends on aspect-ratio and is determined by using 
figure 2. 

NBRR Number of constant x-locations where ACp*s are to be 

interpolated (limited to twenty-five). 

NCOHTS = 0,. Initial locations of free elements will be 

calculated in the program. 

= 1, Initial locations of free elements will be read 
from data cards. 

MITER Maximum number of iterations to be performed 

(usually between 10 and 15) 

IPUNCH = 0, Coordinates of free elements will not be punched 

out after last iteration. 

= 1, Coordinates of free elements will be punched out 
after last iteration. 

Card 3. Format (6F10.5) 

XXL(l) Leading-edge x-coordinate of the root chord. 
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I 


XXT(l) 

YLd) 

XXL(2) 

XXT(2) 

YL(2) 


Tralling-edge x-coordinate of the root chord, 
y-coordinate of the root chord 
Leading-edge x-coordinate of the tip chord. 
Trailing-edge x-coordinate of the tip chord. 


y-coordinate of the tip chord. 


Card U. Ponnat (7F10.5) 


ALPHA 


Angle of attack (in degrees). 


AMACH 


Mach number. 


DELTA 


Length of a segment of leading-edge free vortex 
elements (may be taken as 0.15 C ). 

Length of a segment of wake elements (may be taken 
as 0.15 Cp). 

x-coordinate beyond which free elanents of leading- 
edge- and wake system are represented by a single 
element going to infinity. 

Reference chord. 


Total reference wing area. 


Card 5 . Format ( 8 FIO. 5 ) 


XBRR(I), 


Constant x-locations where 


1=1, NBRR ACp’s are to interpolated. 


Card 6. Format (8F10.5) 


CTT(I), 


Sectional leading-edge thrust coefficients for 


I = 1,(NSW-1) spanwise strips. All these values are set equal to 
zero for complete leading-edge separation. 

*** If NCONTS * 0, go back to card number 1 
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Csrd 7 . Fonaat (1012) 

REyi(l)» (tee higlMr thaa the nuaher of eepimts for each leadlng- 
I ■ 1,(HSW-1) edge free vortex element (numbered from root to tip), 
(terd 8. Format (8F10. ?) 

XE(k) x-coordinates of the end-points of segs»nts of Ith 

K ■ 1,NEIM(I} leading-edge free vdrtex elownt. 

(terd 9* Format (8 f10. ?) 

YE(K), y-coordinates of the end-points of segments of Ith 

K B l,NEUI(l) leading-edge free vortex element. 

Card 10. Format (8F10.5) 

ZE(K) z-coordinates of the end-points of segments of Ith 

K = l,NEli4(l) le8Miing-edge free vortex element. 

•••Cards 8 thru 10 are repeated (NSW-l) times.*** 

Card 11. Format (1012) ’ • 

NNEIM(I), One higher than the number of segments for each wake 
I = 1,NSW element (numbered from root to tip). 

Card 12. Format (8F10.5) 

XXE(K) x-coordinates of the end-points of segments of Ith 

K = 1,NNEIM(I) trailing wake element. 

Card 13. Format (8F10-5) 

YYE(K) y-coordinates of the end-points of segments of Ith 

K 1,NNED4(I) trailing wake element. 

Card lU. Format (8F10.5) 

ZZE(K), z-coordinates of the end-points of segments of Ith 

K = 1,NNEM(I) trailing wake element. 


*** Cards 12 thru lU are repeated NSW times.*** 


**• Go back to card nuaber !.••• 

Note: The punched data cards obtained by running this program with IPUNCH ■ 1, 

ccui be directly used for cards 7 thru for further iterations. 
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k. OUOTUf DATA FORMAf 


All tb« input data osnrds for each case are listed at the beginning of the 
output. The output data at each iteration step is as follows: 

The title card (input data card number l) is printed-out as it is inputted. 
The angle-of-attack (in degrees). Mach Number and iteration number are also 
listed. The end-point locations of the leading-edge free elements are listed 
next. The first row of numbers in each group are the x-coordinates, second 
row tte y-coordinates and third row the z-coordinates. The end-point locations 
of the wake elosents are listed in the similar manner. On the next two pages, 
x-y and y-z digital plots for leading-edge free-elanents are made. It is to 
be noted that the leading-edge elements lying in the plane of the wing and 
along center line are not plotted. So, the elenents next to the center line 
are represented by "l". When similar numbers are connected by straight lines, 
they represent the path of a free vortex element. A "+" sign represents a 
duplicate point. In these two plots there are (NSW-2) rows of free elements. 
The digital plots for wake elements are made on next two pages. The elements 
along center line are again not plotted. There are (NSW-l) elements. 

Some of the intermediate variables are listed under following labels: 

X/C Percent chord location 

2Y/B Percent span location 

GAMAY Bound vortex density over the wing (v ) at the given 

y 

(X/C, 2Y/B) 

CAPGAMA Strength of leading-edge free element (F) at the given 

2Y/B 


DELTA-CP The total ACp at the given (X/C, 2Y/B) 
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The seetional properties are listed under the following labels: 

I Spanwise station nwnber (nuBbercd from root to tip) 

CLI The sectional lift coefficient 

CMI The sectional pitching moment coefficient about the 

y-axis 

CDI The sectional induced drag coefficient 

CTI The sectional leading-edge thrust coefficient. 

The total lift, pitching moment, induced drag and leading-edge thrust 
coefficients are listed after sectional properties. The spanwise pressures 
at constant x-locations are listed under following labels: 

Y y-coordinate 

2Y/B( LOCAL) Percent span location based on local span 
DUiTA-CP The total ACp at the given (x,y) 

The last item listed for each iteration is the absolute force acting on 
leading-edge free elaaents. 

The last page of the output is the "Summary Sheet", which is used to pick 
up final converged solution. It has the following format: 

The title (input data card number l) is printed again. The angle of 
attack (in degrees) and Mach Number are also listed. The other variables 
listed are, 


ITERATION 
CL • 

CM 

CD 

CT 

GMSUM 


Iteration number 

The total lift coefficient 

The total pitching moment coefficient about y-axis 
The total induced drag coefficient 
The total leading-edge thrust coefficient 
Total sum of the strengths of leading-edge free 
vortex elements, except the one at the center line. 


PllSUt P«rcent change in (^eUN values of tvo consecutive 

iterations 

TFA^ fotal absolute force acting on leading-edge free 

elements 

This program has irat yet been coaq>letely automated and the converged 
solution is to be picked hy the user, from the SuBanajry Sheet, by using the 
following criteria: 

The solution is assumed to )mve converged if in two consecutive 
iterations the difference between the total strengths of leading-edge 
free vortex elements is less than one percent and the absolute force acting 
on the free elmnents is in the neighborhood of a minimum. 
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6. APPETOICB 


6.1 APPEiDIX Ai E3MMPLE IHPUT ABD OUTPOT 

foXloviag ii an ucaaple of delta ving of Mpect-ratlo 2 at an angle- 
of-attaek of 30 degrees. The flov is as8\amed to be c^pletely separated 
from the leading«ed^ a»l so sectional leading-edge thrust coefficients are 
set to zei*o on card ntssber 6. Listing of input data cards is ftiven below: 


Listing of Input Data Cards 


6M0 NynMR 


COIUHn Nuntta 

mmiimm2ti2m)isimit4«M4««4««s«isfss»f»*6*666*(>66TmTrT77Ti 

in4y«Tl«Ol21M6TMOlIl«S*7<«OUl«S«7MOtll*t67l«01Z14*«7t99m4t»7«901I}<>}67«90 


iSaiCT •9TI0 • 2.0 
* 7 » 0 10 0 


0. 

9. 

0. 

4. 

9. 

2. 


10. 

1. 

0. 

2. 

0.0 

1. 

0.0 

1.9 

7.9 

1.79 

2. 

0. 

0. 

tND 

0. 

0. 

0. 

0. 

0. 



Output data is listed on the following pages. An inspection of the 
"SuDPnary Sheet" suggests that the converged solution has been reached at 


8th iteration. 
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6.2 APPENDIX Bi COMPUTER PROGRAM LISTING 


A listing of the computer program is given on the following pages. 
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